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The origin of ai^ and e^' orderings in Na2;Co02 

D. Pillay,^ M.D. Johannes,^ I.I. Mazin,^ and O.K. Andersen^ 

^Code 6393, Naval Research Laboratory, Washington, D.C. 20375 
^ Max-Plank-Insitut fiir Feskorperforschung, Stuttgart, Germany 

It has often been suggested that correlation effects suppress the small e'g Fermi surface pockets 
of Nai:Co02 that are predicted by LDA, but absent in ARPES measurements. It appears that 
within the dynamical mean field theory (DMFT) the ARPES can be reproduced only if the on-site 
energy of the e'g complex is lower than that of the Oig complex at the one-electron level, prior to 
the addition of local correlation effects. Current estimates regarding the order of the two orbital 
complexes range from -200 meV to 315 meV in therms of the energy difference. In this work, we 
perform density functional theory calculations of this one-electron splitting A — taig — te' for the full 
two-layer compound, Na2i,Co204, accounting for the effects of Na ordering, interplanar interactions 
and octahedral distortion. We find that taig — te' is negative for all Na fillings and that this is 
primarily due to the strongly positive Coulomb field created by Na"*" ions in the intercalant plane. 
This field disproportionately affects the aig orbital which protrudes farther upward from the Co 
plane than the e'g orbitals. We discuss also the secondary effects of octahedral compression and 
multi-orbital filling on the value of A as a function of Na content. Our results indicate that if the 
e'g pockets are indeed suppressed that can only be due to nonlocal correlation effects beyond the 
standard DMFT. 



Since shortly after the discovery of superconductivity 
in hydrated 1.4H2 O-Nao.asCoOai, there has been con- 
troversy surrounding the discrepancy between the cal- 
culated and observed Fermi surfaces of the parent com- 
pound, Naa;Co02 {x = 1/3). Density functional (DFT) 
calculations show a large central FS, surrounded by six 
small elliptical hole pockets^, while angular-resolved pho- 
toemission (ARPES) measurements find the large pocket 
but do not see the smaller ones^ii^i^. The presence or 
absence of these small pockets drastically changes the 
response properties of the system, and a number of pro- 
posed superconductivity models depend on theni^i^i^ii^. 
There are two basic ways to resolve the controversy: ei- 
ther the experiment is not representative of the bulk 
electronic states, or correlation effects beyond the DFT 
calculations qualitatively change the Fermiology. The 
latter seems rather plausible, given the narrow width 
of the Co t2g bands and the fact that DFT fails to 
reproduce the nonmagnetic ground state for x < 0.7. 
The effects of adding correlation have been recently ad- 
dressed through DMFT calculations which account for 
both on-site Coulomb repulsion (Hubbard U) and local 
fluctuation effects (albeit not for long-wave-length ferro- 
magnetic calculations, compatible with the observed at 
X > 0.5 magnetic interactions in this system). Despite 
initial controversyiiii^ii^, it has been firmly established 
that DMFT calculations only agree with ARPES if the 
centroids of the the aig and e' Co d— bands are such 
that eai is higher than e^' , which would yield a posi- 
tive splitting, A = Eai — Ee' > 0. If, on the other hand 
Eaig — Ee' < then the small hole pockets are preserved 
and even slightly enlarged. This result does not depend 
on the value of U or the particular Hamiltonian used in 
the calculation^^. The sign of A is therefore of great 
importance, but difficult or impossible to access exper- 
imentally. A tight binding fit to full potential density 
functional theory calculations^ employing all 5 Co d or- 



bitals and 3 O p orbitals in a double layer system found a 
strongly negative A ~ -200 meV for x =0.3. Other tight 
binding models fit to a one-layer system and employing 
a smaller basis set of Co t2g orbitals^^'^*^ found A ^ -144 
meV, and an alternate procedure for extracting A from 
the band structure yielded a value of ~ 100 meV^^. On 
the other hand, quantum chemical calculations of em- 
bedded clustersi^iii invariably find a positive A = 315 
meV. There has so far been little understanding of what 
might cause such strong discrepancy between the values 
obtained through the two different methodologies. 



In this work, we present DFT calculations of A that 
include both layers of the double-layer unit cell and real 
Na ions (no virtual crystal or other approximation) and 
we extract our values directly from the DFT calculation 
without employing a tight-binding fit. It is crucial to 
understand that A is controlled entirely by the electro- 
static field at the Co site and Co-O overlap integrals, 
which arc both accurately described by DFT. Note that, 
within DFT, the Fermi surface itself hardly depends on 
the sign of A, as long as |A| is moderate, but the effect 
of the local correlation effects (as evaluated by DMFT) 
is drastically different. At every Na doping level and for 
every reasonable value of octahedral compression, we find 
that A is negative. We find that the main cause of the 
ordering between aig and e' orbitals is the positive Na 
potential in the interstitial plane that drives the out-of- 
plane aig orbitals downward in energy, a factor that has 
been overlooked in all previous studies^^. Interestingly, 
we find that octahedral distortion as an isolated factor 
increases A , a reversal of crystal field predictions for 
CoOg octahedra. This highlights the strong role of Na 
to re-establish a negative A and emphasizes that Co-O 
hybridization depends on interactions beyond the CoOe 
octahedra. 




b) 




FIG. 1; (color online) Charge density plots for a) the aig and 
b) the e g' complexes of NaCo02. The plots employ squared 
wavefunctions, produced using full DFT calculations, taken 
at the r point where the aig and Cg' states are pure. The 
sum over both squared eg' orbitals is shown in b). 



A. Calculational Methods 

All calculations have been performed using a plane 
wave density functional theory code, the Vienna Ab- 
initio Simulation Package (VASP)i^. The projected aug- 
mented wave method (PAW)-° was used along with the 
generalized gradient approximation^^ to the exchange 
correlation potential. A \/3 x \/3 cell that preserves 
the overall hexagonal symmetry was used to simulate the 
bulk. Ground state geometries were fully optimized us- 
ing a Monkhorst Pack k-point mesh of 4x4x2. For the 
total energy calculation and density of states (DOS) the 
Brilliouin zone integrations were performed on a 8x8x4 
mesh. 

To calculate A, we first consider the Co plane in a 
cubic coordiate system where three Co atoms are located 
equidistant along the x,y, and z axes, i.e. they form a 
plane perpendicular to the [111] cubic axis which is also 
the z axis in the hexagonal system. In this cubic system, 
it is easy to explicitly write out the aig and e' orbitals: 



XY + YZ + XZ 



Hg 



XY 



XZ -YZ 



(1) 



V3(YZ- XZ) 



Standard band structure codes project the density of 
states onto the three symmetry representations appropri- 
ate for a hexagonal symmetry, that is, Nz2 projects onto 
the 3z^ — 1 orbital, and N^2_y2^^y and Nxz,yz, each give 
the total projection onto the two orbitals belonging to 
each respective representation. Since the projected den- 
sity of states contains the square of the projected orbital, 
we square each term in Eq. [U keeping only the terms con- 
forming to the correct symmetry. The resulting recipe is: 



Na,^{e) = N^e) (2) 

A^e',(e) = N,x2_y2^xy{t) + 2Nxz,yz{e) (3) 

From these expressions, we calculate band centers of 
each complex. It is important to note that the e' states 
share the same irreducible representation as the unoccu- 
pied Cg states and that there is mixing between the two 
complexesi^i. Therefore, without specifically projecting 
out only the e' states as in our methodology, erroneous 
values for A will be found. The onsite parameters that 
result from this procedure are not identical to the onsite 
parameters in the effective 3x3 Hamiltonian, so far used 
in DMFT calculations because the latter contain not only 
the Co-based aig and e'g states, but also the hybridized O 
states. The latter can be obtained only by downfolding 
the calculated Bloch functions onto Wannier functions, 
as has been done befor o^^i^^i^^i^^ . Still, our method can 
be used to estimate the changes in A brought about by 
different factors, such as octahedral distortion, Na fill- 
ing, or Na ordering. In the several cases where we have 
compared our values of A to those calculated from the 
Wannier functions of a three band model, we have found 
that the latter are always more negative than the val- 
ues extracted directly from DFT. It is worth noting that 
these Wannier functions are very long range, having a 
sizeable weight not only on O and on nearest neighbor 
Co's, but even at the third neighbor Co sites. This may 
have an effect on the onsite energy ordering, although 
we believe it will be quantitative, not qualitative. An 
additional effect of this extent is a strong reduction of 
the Hubbard U as compared to typical values for local 
d-orbitals in transition metal oxides. 



B. Results and Discussion 

We first address the effect of the octahedral compres- 
sion (trigonal distortion) that lowers the symmetry at 
the Co site and splits the t2g triplet into the aig and 
e' complexes. Some previous works have assumed that 
the distortion destabilizes the aig state^^i^ii^i^^, while 
others assume oppositely that the e' state is higher in 
energy^2i21. We have examined this question from two 
points of view. First, we have constructed a tight-binding 
model of a CoOg octahedron, including the three Co 
t2g orbitals and three p orbitals on each of the six O 
atoms. We have diagonalized the resulting matrix for 
both compressive and expansive distortions. Using the 
canonical scaling of the tdpa and tdp-n parameters, namely 
idpa = yitdp-m we find that the compressive distortion, 
as seen in Naa;Co02 for all x, destabilizes the e' state. 
However, if the ratio of the two parameters were to be 
increased by approximately a factor of two, as has been 
seen in some casea^^, the order of energies is reversed and 
the aig state is destabilized. Second, we have performed 
a full DFT calculation of A as a function of the trigonal 
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FIG. 2: (color online) A = ea^^ — e^i as a function of octa- 
hedral distortion in the jellium filled and real Na filled com- 
pounds. The vertical bar at ^ 5.2% marks the position of the 
octahedral compression obtained through relaxing the oxygen 
height in our calculations. 



FIG. 3: (color online) A = ea^^ — te' as a function of electron 
filling for jellium filled and real Na filled compounds. The 
critical doping Xc where the eg' complex becomes completely 
filled is marked by a vertical bar. 



distortion. We use a fictitious compound, e"Co02, in 
which one electron per formula unit has been added to 
the C0O2 compound to achieve the band-insulating state, 
analogous to NaCo02. The extra charge is compensated 
by a uniform positive background. We use this technique, 
termed "jellium filling" , to eliminate the effects of the Na 
potential and of unequal band fillings, which we discuss 
in later section. As can be seen in Fig. [51 octahedral 
compression always produces a positive A, i.e. the aig 
state is raised above the e' . This demonstrates that the 
scaling of t^pa and tdp-n is outside the canonical range and 
points to the importance of neighbors more distant than 
the surrounding O atoms. 

By far the strongest factor determining the value of 
A is the presence of ionized Na atoms in the intercalant 
plane. The positive Coulomb field they produce serves 
to strongly shift the aig orbital downward compared to 
the e' orbital. As seen in Fig. [U the aig orbital has a 
greater extent along the z direction (c-axis of the crystal) 
and therefore penetrates further into the Na plane. The 
resultant lowering of the aig states is sufficient to shift the 
sign of A from positive to negative. Replacing the jellium 
filling discussed above with real Na ions (again one per 
formula unit), we find a negative shift of approximately 
150 meV in A for all trigonal distortions, as seen in Fig. 
[51 Because quantum chemical calculations account only 
for the filling effects of Na and neglect the electrostatic 
field produced by real ions, this shift is not present and 
A is consequently found to be strongly positive. Not 
surprisingly, when Na ions are located directly atop a Co 
ion (26 site), the effect is stronger than when they are 
located equidistant between three Co ions (2c site). 

The final factor contributing to the size of A is the 
non-uniform filling of the two band complexes. For all 
levels of Na intercalation, the aig bands are less filled 



than the e' bands. At a critical filling of Xc 



0.67, 



the e' complex becomes entirely filled, leaving holes only 
in the aig bands. The increased on-site Coulomb repul- 
sion due to greater filling drives the e' upward in energy 
compared to aig bands. In the absence of real Na ions 
(hence no positive interstitial Coulomb field) and with- 
out octahedral distortion, the calculated A for C0O2 is 
—77 meV, refiecting this multi-band effect. We have cal- 
culated A for Naa;Co02 and for a;e~Co02 (jellium fill- 
ing) as a function of x, with all oxygen positions fixed 
to perfect octahedra. For the jellium filled compounds, 
A initially remains approximately constant and then in- 
creases dramatically at Xc '^ 0.67 where the e' bands 
become entirely filled, as seen in Fig. [31 The Na:rCo02 
curve, on the other hand, decreases initially and then 
fiattens out once the e' complex is full. In the former 
case, electrons are added to both e' and aig states at low 
X, pushing both complexes upward nearly equally. Once 
the e'g holes are filled, all further electrons enter the aig 
complex, pushing the aig onsite energy up relative to e'g 
and causing A to increase. The situation is more com- 
plex when real Na ions are added. At low x, both e' and 
aig, electrons are added to both complexes and on-site 
repulsion raises both energies together, but the increas- 
ingly positive Coulomb field created by adding Na ions 
pushes the aig states downward. A therefore increases 
with X as long as holes remain in both symmetry com- 
plexes. After Xc, where additional electrons enter only 
the aig complex, the positive Coulomb field acts oppo- 
sitely to the increased on-site repulsion. Apparently the 
magnitude of these two effects is approximately equal, as 
evidenced by the nearly constant A at a; > x^,. 

The effect of octahedral distortion on A competes with 
the effect of the positive Na"*" potential, and both are a 
function of the Na filling, x, which separately affects the 
aig and e' ordering. Having studied each of these fac- 



TABLE I: A = eaig - te' values for full DFT calculations, 
including relaxation of the lattice parameters and ionic posi- 
tions. Values for Na placed at Nal {2b trigonal site) and at 
Na2 (2c octahedral site) positions are given for comparison. 

X Nal position Na2 position 






-80 meV 


- 


1/3 


-100 meV 


-95 meV 


2/3 


-114 meV 


-101 meV 


1 


-77 meV 


-6 meV 



tors individually, we proceed to calculate A as a function 
of X, accounting for all factors simultaneously. This is 
accomplished by relaxing the lattice parameters (both a 
and c) and the ionic positions at each Na filling and us- 
ing supercell calculations to include real Na ions in the 
interstitial plane. Our results are tabulated in Table [J 
Lechermann et al have derived A for x=0.5 with Na in 
the Na2 position, by integrating the partial density of 
states and also by constructing first principles Wannier 
functions, arriving at an estimate of -100 meV, in ex- 
cellent agreement with our Table 1. Their estimate for 
Delta at x=l is also consistent with oursi^. It is clear 
that A is a non-monotonic function of x, owing to com- 
petition between the various effects associated with Na 
filling, but that it remains negative for all Na content 
< a; < 1. 



C. Conclusions 

The crystal field splitting A tai^ — f-e' is known to be 
the parameter that controls the presence or absence of the 



controversial e' pockets at the Fermi level in Naa;Co02. 
Here we show that, within the DFT framework which 
is highly reliable for calculating one-electron parameters, 
A is always negative and strongly so for low x. We have 
shown that this fact is mainly due to the electrostatic 
crystal field of Na as well as demonstrating that the O 
ligand field (octahedral distortion), in fact, leads to the 
opposite sign of A. The explicit effect of ionized Na has 
been neglected or severely approximated in previous cal- 
culations, leading erroneously to positive values of A. 
This result indicates that local correlation and/or fluc- 
tuation effects are not responsible for suppressing the e' 
Fermi sheets. Either longe-range fluctuations, outside 
the scope of DMFT, or strong deviations from the bulk 
structure/stoichiometry such as surfaces, defects, or va- 
cancies, must therefore account for the non-observation 
of these pockets in ARPES experiments. 
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